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Abstract 

Phenomenological analysis to the R-parity violating supersymmetry with a vector-like extra 
generation is performed in detail. It is found that, via the trilinear couplings, the correct neutrino 
spectrum can be obtained. The Higgs mass rises to 125 GeV by new up-type Yukawa couplings of 
vector- like quarks with no need of very heavy superpartners. Phenomena of new heavy fermions 
at LHC are predicted. 
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I. INTRODUCTION 



Recently a standard model (SM) Higgs-like particle with a mass of 125 — 126 GeV was 
discovered [l|. In the paradigm of the weak scale supersymmetry (SUSY) which aims at 
the naturalness of the electro-weak scale, however, such a Higgs mass brings in tensions, 
especially the minimal SUSY SM (MSSM). Nonminimal and still natural scenarios of SUSY 
are thus motivated. One of them is the MSSM with a vector-like generation gives 
the right Higgs mass naturally, is consistent with precision electroweak measurements, and 
has a rich phenomenology {2-6]. In the framework of SUSY, vector-like fermions can also 
be motivated by other theories beyond SM, such as SUSY extension with extra-dimensions 
or with composite states [7]. So it is worth asking the question whether such a scenario also 
provides explanations to other problems such as neutrino masses. ^ 

Neutrino oscillations are the undoubted new physics beyond the SM. Daya Bay |^] and 
RENO jo] experiments recently discovered a relatively large ^13 ^ 8.8° ± 0.8°. Within the 
framework of SUSY, in the absence of R-parity conservation, neutrino masses and mixings 
can be generated from lepton number violating (LPV) couplings 10|. This approach was 
extensively studied before ll| . It is known that all the neutrino experimental results, includ- 
ing that of oscillation phenomena like the large atmospheric mixing angle ^23, the hierarchy 
of oscillation frequencies Am^i ^ ^77132 and the smallness of ^^13, can be understood in 
three generation LPV MSSM. However, this needs some special requirements for relevant 
coupling constants and mass parameters. 



Combining bot 
extra generation 



1 considerations above, we will work in the LPV MSSM with a vector-like 
^. While this model takes the vector-like slepton doublets as the two 
Higgs doublets needed for the electroweak symmetry breaking, the SM-like Higgs mass can 
be naturaly 125 GeV ^. Extra trilinear LPV couplings between ordinary fermions and 
vector-like fermions provide a much larger parameter space to explain neutrino pheomena 
right. 

In this paper, phenomenological aspects of the model will be analyzed. In Sect. II, we 
make a brief review of the model. In Sect. HI, neutrino masses are calculated. For the 
neutrino physics, noting the enlarged parameter space, we consider trilinear LPV couplings 
carefully. One-loop contribution to neutrino masses due to new trilinear LPV couplings 
is calculated, theoretical analysis are performed and numerical results are shown in detail. 
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Besides, we analyze the SM-like Higgs mass and explicitly show that it can be increased to 
125 GeV by two new Yukawa couplings of the up-type Higgs with vector-like quarks in Sect. 
IV. The LHC phenomenology of the new fermions is analyzed in Sect. VI. The summary 
and discussions are given in the last section. 



II. A BRIEF REVIEW OF THE MODEL 

This modelj4| is SUSY and SM gauge invariant, and R-parity violation with baryon 
number conservation is assumed. For the matter content, in addition to the ordinary 3 
generations (3G), a vector-like generation is introduced in. Without R-parity conservation, 
this can be also thought as that there are 4 + 1 chiral generations, where '4' stands for four 
chiral generations with SM quantum numbers and '1' for another chiral generation with 
opposite quantum numbers. The 4 chiral generations with same quantum numbers mix. 
The '1' has Dirac masses with only one combination of the '4', thus, there are always SM 
required three massless chiral generations and one massive vector-like generation. 

In terms of mass eigenstates (before electroweak symmetry breaking), the massive slep- 
tons in the vector-like generation are taken as the two Higgs doublets. New particles beyond 
the MSSM are the following with quantum numbers under SU(3)cXSU(2)ixU(l)y, 

^4^(1, 1,2) , E|,(l,l,-2) , g4(3,2, I) , QH(3,2,-i) , 

f/4^(3,l,-|) , f/]^(3,l,|) , D2{3,1,1) , D^(3,l,-|) . 
The superpotential is conveniently written as 

>V = Wo + , (1) 

where Wo and stand for that with lepton number conservation and LPV, respectively. 



+yfQ^H^Uf + y'^'^Q^HuUl + y^'QuHdUfj + y'^'QuHuD^ , 



and 



W^ D XijkLiLjEl + X'^jf.QiLjD'j. + XfjLiLjE^ + XfjQ^LiD'j 

+XgQiLjDl + Xf^'Q^LiDl + X^QhUU^j . (2) 
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where Lj, Qj, Ef, D^, Uf, i=l-3, are the first three generation SU{2)l doublet leptons, 
doublet quarks, singlet charged leptons, singlet down-type quarks and singlet up-type quarks, 
respectively. and are the up-type and down-type Higgs. Note that the term Q^HuD^ 
in Wo was missed in Ref. j4|.^ And in interactions of purely singlets are omitted, which 
are irrelevant to our study. 

By assuming universality of the mass-squared terms, the alignment of the B terms the 
soft mass terms and the trilinear soft terms of all fermion's superpartners in the model are 

-CD M^L\U + M^H\H, + MlHiHu + MIE^E'}^ + M^Q^nQm + M^[/;;t;7c^ 

+Mlb'iD'i^ + Ml^E^-E^^ + MI^Q^^Qh + M^^U'SU'h + Ml^D'^D'k (3) 
+ {BfiHdH^ + B'u'EIEIj + B'^fi'iQ^QH + B"" ^i^'UIU'h + B^'^i^'DID'h + h.c.) . 

Proper values of the new B'^'^'^ ^'^'^'^ terms are set to avoid unwanted color symmetry 



and purely U{1)y symmetry breaking, see Eq. (11, 12) in paper {4], therefore EWSB in our 
model is just the same as in MSSM. After EWSB, the specific fermion mass matrixes and 
sfermion mass-squared matrixes are given in Appendix A. 



III. NEUTRINO MASSES AND MIXINGS 



LPV results in nonvanishing neutrino masses. In this model, in addition to traditional R- 
parity violation in the MSSM, a lot more bilinear and trilinear LPV interactions are brought 
in through the vector-like generation. In this work, the trilinear R-parity violating interac- 
tions will be studied. To avoid complication due to too many LPV sources, sneutrino VEVs 
will not be considered. There are several reasons for this. First, we can phenomenologically 
assume the universality of the soft SUSY breaking mass terms at the weak scale, to avoid 
dangerously large flavor changing neutral currents (FCNCs), without considering any UV 
completion of the model. In that case, because of the alighnment in bilinear terms of the 
superpotential and that of soft terms, R-parity violating bilinear terms can be rotated away 
via field redefinition, and sneutrino vacuum expectation values (VEVs) vanish in the phys- 
ical basis. The second reason is from consideration of underlying models. SUSY breaking 

^ It modifies the down-type fermion mass matrix and scalar mass-squared matrix. Correct ones, as well as 
the resulting mixing matrix are given in the Appendix A. 
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is introduced effectively in our model, it can result from gauge mediated SUSY breaking. 
Then the messenger scale can be as low as 100 TeV, even if the universality scale is at the 
SUSY breaking messenger scale, the running effect is small, and the bilinear LPV is not 
important compared to the trilinear ones. Finally, small sneutrino VEVs can be included 
in the analysis nevertheless in future works, after the role of new trilinear LPV interactions 
gets a thorough understanding. 

The trilinear LPV Lagrangian relevant to neutrino masses is from Wfi, 

C C -Xijkil*kRi^iRljL + IjihR^iL) - Kjki^lRi^iRdjL + djidkR^iL) 

-Xf^iElutRhL + hiEfv^L) - X^idlRiy^RQA + QJkRi^iL) (4) 
-\f;{Dlu^Rd,L + d,LDfu,L) - Xf^'iDMRQA + Q^Dfu^L) 
-X^iUn'^-RQH + QuUiyrL) + h.c. . 

where v^j^ stands for the left-hand neutrino. 

The 7 types of trilinear LPV interactions in the above equation induce 14 types of one-loop 
diagrams contributing to the neutrino spectrum, which are proportional to AA, A'A', A^A^, 
AA^, A«A«, A«A^, A^A^, A^A^, A«^A«^, A'A^, A'A^, A'A«^, A^A^^, A^A«^, respec- 
tively. The Feynman diagrams and the corresponding analytical results are shown in Fig. 1 
in Appendix B. For simplicity and without losing our purpose, in the Yukawa interactions of 
Wo we assume that only , y^' , y'^' , y^, y^ , y^ , y^' are nonvanishing, that is vector-like 
particles have Yukawa interactions only with the third generation. Thus, the vector-like 
generation has little constraints from the collider phenomenology. 

Before starting to analyze the neutrino mass spectrum, some assumptions are intro- 
duced in order to control the parameter space and get relatively simple analytical result. 
Since four new up-type Higgs Yukawa couplings y^ , y^, y^'^, y^' and five new down- 
type Higgs Yukawa couplings y^, y^, y'^' , y^^, appear in our model, and among 
which y'^^, u'^^ ■, ■, U^' provide the mass mixings between vector-like generations, and 
further more, they have infrared quasi- fixed point [sl, so we assume i/*^^ = y^ = and 
yQU ^ yH' = y*, < 1. We also set = y^' = 0, y^ < 0.04 and y^ ^ = < 0.08. In 
other words, we neglect all new down-type Higgs Yukawa couplings in quark sectors while 
consider all of the new up-type Higgs Yukawa couplings only and take <^ which is a 
reasonable assumption. 

Basing on the above assumptions, contributions from AA, A'A' type diagrams can be 
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simplified to the familiar forms 12|-|l4| 



1 Tr 

^^^Ua - ^Ai33Aj33m^sina^ cosof In— , 

3 & b 

Mijlx'y ^ -^{Kss^'jss^b sin cos In ^ + A-asA^^am, sin cos In 3 

+ A^32A'-23 sin ag cosag In (5) 

sl 

where in the first equation, we only keep the dominant contributions and in the second 
equation, we keep the dominant and subdominant ones, a,-, a^, as, at are the angles of 
the corresponding 2x2 fL{ji),bL(^R), SL{R),tL(R) unitary matrices. Unfortunately, the other 
equations, (B3)-(B14) in Appendix B, can not be simplified by following similar process 
because there are mixings between different vector-like generations. So these can only be 
analyzed numerically and will be discussed later. 

At last, without loss of generality, among all 7 types of LPV trilinear couplings, we take 
4 of them, A^, A^, A*^ and A^, for consideration while assuming the rest of them. A, A' and 
are negligible. The realization through different LPV trilinear coupling combinations 
can be derived straightforwardly. The method to calculate the neutrino mass matrix we use 
is given in Appendix C. 

Here we list the parameters of neutrino oscillation given by experiments, Arn^i = (7.59 ± 
0.21) X 10-^ eV^ Aml^ = (2.43 ±0.13) x 10"^ eV^ and sin^ 2^12 = OMlto^oH sin^ 2^23 > 
0.92, sin^ 2^13 = 0.088 ± 0.008. Scanning the parameter space with proper EWSB, we find 
by adjusting the ratios and values of the LPV trilinear couplings we choosing, the correct 
neutrino spectrum can be generated through the A^A^, X^X^ and X^X^ type one-loop 
diagrams. Numerical illustration is shown in Table I, in Set I we take the mass mixings 
assumptions mentioned before, in Set II we take different mass mixings and bigger vector- 
like masses for comparison. The specific parameters settings see Appendix C. 

That is by choosing (for Set I) 

13 n oc; 33 A \Q o ^ in-6 \Q \D \Q \D \Q \D 

Tq ~ ^-^^ ' To" ' ^2 ^ ' ^ ^2 ' ^3 ^ ^2 ' ^ ' 

A23 A23 



\ ~ 1.4 , \ ~ 1 



\H - ' - - ' -^2 ~ -^23 > -^13 ~ -^23 ~ "^33 ~ -^23 ' (6) 

^^2 ^^2 



we have 



^ ~ 0.17 , m^2 ~ 5.9 X 10"^ eV , m^i ~ 5.1 x 10"^ eV , 

sin ^13 ~ 0.143 , sin ^23 ~ 0.581 , sin ^12 ~ 0.559 . (7) 







\.3^j3 


'■J '^i3^j3 




^I^-Ia^a^ 


Set I 


0.0043 


0.238 








0.08 


Set II 


0.0027 


0.168 


0.004 


0.003 


0.011 



TABLE I: Numerical illustration for 5 types of one-loop contributions in our model ,tlie specific 
parameter settings see Appendix B. M^- (GeV) stands for the parts in Eq. (4,5) excepting the 
LPV trilinear coupling constants. 

Unlike in the 3G LPV case, where A'j33Aj33, \'i2i,'^'i2,2 + ^^32^^23 ^'^^ Aj33Aj33 type one-loop 
contributions are dominant, subdominant and next-to-sub dominant, here in our model, un- 
der the assumptions mentioned before, A^Aj^, ^f^f and A^A^ type one-loop contributions 
are dominant, subdominant and next-to-sub dominant, respectively. This is because the new 
fermions ri, ti^2, &i,2 in the internal lines, see Fig. 1, are much heavier than the third 
generation fermions r, t, h. 



\ A 



tA 



1,2 



13 



Af / 



tA 



1,2 



H 



FIG. 1: New one-loop contributions to the the neutrino masses and mixings from A^A^, X^X^ 
and X^ X^ type couplings. All particles stay in mass eigenstates. 



For the same reason, our requirements of the new LPV couplings we choose are of order 
10^^ and small enough to avoid measurable FCNC decays such as /i — )■ 67 15|. It worth 
to note that by decoupling the vector-like generation, correct neutrino masses and mixings 
cannot be obtained via AA, A' A' type one- loop contributions. 

In addition, X^X^ type contribution containing up-type (s)quarks in the internal lines is 
absent in 3G LPV models because the vector-like down-type doublet quark mixe with 
the right hand singlet top quark. 

From Table I, we can also see that by choosing A^Aj^, A^A^ and X^Xf^ type one-loop 
contributions, the correct neutrino spectrum can also been generated in parameters Set II, 
we don't list the detailed results here. 
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IV. HIGGS MASS 



There are four new up-type Higgs Yukawa couplings in our model, , y'^, y'^^ , y^' , 
corresponding to the Yukawa mass, ml^^, m\^, m\^, m\j, separately, and also five new 
down-type Higgs Yukawa couplings, y^, y^, y^' , y^^, y^ , corresponding to the Yukawa 
mass, m34, 77134, 77743, '"^44; "^H; separately. The related superpotential contributing to the 
lightest scalar Higgs mass is shown in Wq. According to the assumptions mentioned in the 
last section, we neglect the down-type Higgs Yukawa contributions and the small up-type 
contributions between the SM third generations and the extra vector-like generations. The 
relevant superpotential can be simplified as 

W = fi'^Q^Qn + f^^'DlD'j, + y'^'QuHuD^u + v'^'^QaHuUI , (8) 

So when neglecting the small D-term and the two-loop contribution, the new one-loop 
contribution to the lightest scalar Higgs square-mass is j^, Q 

a™; ^ i^i 1,, _ ^(5 - - i) + 2 ^(1 - J-,1, (9) 

where v = 174Gev indicates the Higgs VEV and 

y^' = = , X = Ml/M'y , V = log^ , (10) 



(Ah' - jJ^H' cot = {Aqu - jJ^Qu cot = {Ay - jj^v cot = X 



in which, for simplicity, /x*^ = = My stands for the vector-like mass of the new up-type 



quarks, Mq = M|, = m'^ (see Eq.(4)) and Ms = \/ My -|- 777^ stands for average mass of the 
new up-type squarks. 

In MSSM, the Higgs mass from the t, t one-loop contributions is about 110 GeV, for 
A* = /i = 400 GeV, mi = 400 GeV and tan/3 = 10. Direct search bounds from CMS 
for exotic heavy top-like quark set limits of M^/ > 557 GeV if B{t' — )■ Wb) = 1 jl3] and 



Mf > 475 GeV if B{t' — )■ Zt) = 1 18|. When considering the mass mixing between the 



vector-like quarks and the SM third generation quarks, in other words, considering the 



realistic branch ratios, the mass limit is adjust to be Mf > 415 GeV [19|, l20|]. So if we set 
the vector-like fermion masses in our model to be My ~ 500 GeV, the soft supersymmetry- 
breaking parameters to be 777 ~ 700 GeV, Ay = iiy ^ 500GeV and Byfiy ~ 500^ GeV^, then 
from Eq. (10), in order to get approximately 125 GeV Higgs mass, for about My = 500 GeV 
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and Ms = 850 GeV, we just need to set yy ~ 1, or say, need to set 



m 



H 



ml 



174 



GeV. These values are just near their infrared quasi-fix point, as mentioned in last section. 
Evoked by the ATLAS and CMS discovery of the enhancement in 77 channel and little 



deviation in ZZ channel 



22 



23j, the effects of the exotic vector-like quarks to the Higgs 



production and decay have been extensively studied recently 



2l| . In general, in a theory 



with vector-li 



21 



ie generations extension, the new fermion contributions are suppressed by 



2J]. So only the very large couplin gs t o the Higgs can obviously enhance 



the Higgs production and decay in the 77 channel 2l|, but as we have mentioned, these 
couplings have quasi-fix point which limits their TeV values to be about 1 5| . This value is 
large enough to accommodate nih ~ 125 GeV, but too small to influence the Higgs decay, 
one can't depend on vector-like fermions by themselves to modify the Higgs decay branching 
ratios. As far as the Higgs problem to be concerned, extra vector-like fermions are mainly 
introduced to adjust the Higgs mass. However, the 77 and ZZ channel anomaly, if they 
persist, can be realized through the light stop scenario [25|, which beyond our scope in this 
paper. 



V. THE EXTRA VECTOR-LIKE FERMION DECAYS 

To be clear, we list the new extra vector-like fermions below: 

in which E"^ mixes with r^; mixes with r/j; Q\,Djj mixes with bi] Dl,Q\j mixes with 
bn, Q\]U'fj mixes with ti, U2,Qh mixes with t/j. These exotic heavy fermions can decay 
into SM bosons, see Fig. 2, which will analyze bellow. Our analysis agree with the results 
given in [5|]- However the slightly difference comes from their neglect of the contributions 
proportional to in the vertex of Feynman rules. 

Note that theoretically speaking, when kinematically allowed, the exotic fermions pre- 
dicted in our model have the other two decay modes: through supersymmetric gauge kinetic 
interactions or the supersymmetric Yukawa interactions, decay into chargino/neutralino and 
sfermions, such as ri — )■ C+tV; bi — )■ Nib, ti — t- C^b , where Ni, i=l-4, is neutralino, is 
chargino; through LPV interactions, see Eq. (2), decay into fermions and sfermions, such as 
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FIG. 2: Tree-level decay of new exotic fermions in our model, all fermions stay in mass eigenstates. 

Tl — )■ e/i, bi —7- ft, tl — i- eb. Although the kinematical conditions for the latter decay mode 
are easy to be satisfied, but we have already seen in section III, the LPV couplings in our 
model, in order to explain the neutrino spectrum, are of order 10~^, so we can neglect this 
kind of decay channels reasonably. On the other hand, for simplicity here in our work, we 
assume the former decay mode is not kinematically allowed. Therefore, the exotic fermions 
can only decay into SM bosons. 



A. Tl decays 

The weak bosons interaction Lagrangian to r, ti is 

C D gV^^uJiL^^rLW- + gl^^jiLi'r^Z^ + gl^^^^fiR^^TRZ^ 

+gtr,fi,Tnh' + CtlTik/i^ + h.c, (12) 

the couplings and the decay widths of ti are given in Appendix D. 



The main characteristic of the lepton sector is that there must be mass mixing between 
the third and the vector-like lepton, otherwise the new heavy charged leptons ti will be 
stable and give unacceptable cosmological heavy relic 26|. For Specific, when = 0, the 
off-diagonal elements of L'^, are equal to zero. That's why we set ^ in section 
II while discussing neutrino spectrum, more specifically, we set < 0.04. Under these 
parameters settings, numerical results of ti decay into W, Z, are shown in Fig. 3, we 
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800 



MtI (GeV) MtI (GeV) 

FIG. 3: The decay widths of the new lepton ri (left panel) and its branching ratios (right panel) 
with = 0.04. 

can see in the limit of large the branching rations are BR(ri -^Wur) ~ 0.7 and 
BR(ri Zr) = BR(ri -» hr) ~ 0.15 . 



B. ti 2 decays 

The weak bosons interaction Lagrangian to t, ti, t2 is 

^ ^ attMYb^W- + g^.MYbLW^ + g^^^JiR^bRW;; + g^^,j2RYbnW;; 
gZbj2Ll'b,LW- + g^^,j2Rl'hnW- + g^^J^Ll^^Z, + g^^J^Ll'tLZ, 
+gi,tj2LYtiLZ, + g?^,J,nYtRZ, + g?^j2RYtRZ, + g?^,j2RYtiRZ^ 
+gtjiLtRh' + gl.jLtmh' + gttMRh' + gttjLt2Rh' 
+9ttj2Ltmh' + gf^tJiLt2Rh' + h.c, (13) 

the couplings and the decay widths of ti^2 are given in Appendix D. 

As mentioned in section II, we take ^ < 0.08, y^^ ~ y^' < 1, the numerical 
reasults are shown in Fig 4, 5. We can see in the limit of large Mt^j, the branching rations 
of ti are BR(ti Wb) ~ 0.4 and BR(ti Zt) = BR(ti h°t) ~ 0.3, the branching 
rations of t2 are BR(t2 ^ Wbi) ~ 0.85 and BR(t2 ^ Zti) ~ 0.15 . 
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FIG. 4: The decay widths of the hghtest new up- type quark ti (left panel) and its branching ratios 
(right panel) with y^^ = = = 0, = 0.08 and y^^ ~ y^' = 1. 
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Zti 



FIG. 5: The decay widths of the heaviest new up-type quark t2 (left panel) and its branching ratios 
(right panel) with y*^^ = y^ = yD = yU ^ yQ ^ 0.08 and y^'^ ~ y^' = 0.98. 
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C. 612 decays 



The weak bosons interaction Lagrangian to b, hi, 62 is 

+9i^bj2LYhLZ, + gl^.linYbRZ^ + gl^,l2RYbRZ^ + gl^,j2Rl^biRZ, 



+C6i«^2L&li?/^^ + gi^,JlLb2Rh'' + h.C. 



(14) 



the couphngs and the decay widths of bi^2 are given in Appendix D. 




300 400 



500 600 700 
Mb2 (GeV) 



800 900 1000 



300 
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FIG. 6: The decay widths of the heaviest new down- type quark b2 (left panel) and its branching 
ratios (right panel) with y^^ = = = ^ yQ = 0.O8 and y'^^ ~ y^' = 0.98. 



The numerical results of b2 decay widths and branching rations are shown in Fig. 6 
under the parameter settings mentioned before, we can see BR(b2 — )■ Wbi) ~ 0.85 and 
BR(b2 — )■ Zbi) ~ 0.15. The branching rations of bi are BR(bi — )■ Wt) = 1 which are not 
shown here. 
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VI. SUMMARY AND DISCUSSION 



We have studied several phenomenological aspects of the LPV MSSM model with a vector- 
like extra generation: the neutrino spectrum, the Higgs mass and the LHC phenomenology 
of the new predicted fermions. The results are: 

• The correct neutrino masses and mixings, especially the relatively large ^13 can be 
generated from trilinear LPV couplings. The new trilinear R-parity violating couplings 
make it easy to generate the proper value of ^13. These couphng constants need to be 
about 10"^. 

• The two new up-type Higgs Yukawa couphngs, y^' and y^^ ., between the vector-hke 
quarks and the SM third generation quarks, with values about 1 near to their infrared 

quasi-fixed point in TeV scale, can give rise to 125 GeV Higgs mass with no need of 
very heavy new superpartner. 

• There are five new heavy fermions, ri, ti 2, bi 2 , predicted in this model. They can only 

decay into SM bosons by some kinematic assumptions. The branching radio depend 
on the mass mixing between the vector-like fermions and the SM third generation 
fermions. These charged exotic fermions would be quasi-stable if such mass mixings 
are very small. 

Based on our previous work about bilinear LPV couplings, further research on the renor- 
malization group (RG) of them is worthy to be studied in the future. There are also plenty 
of aspects to be further analyzed in the area of new fermion LHC phenomenology based on 
this model. 
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Appendix A: THE (S)FERMION MASS MIXING MATRIXES 



Because the CP violation is not considered in this paper, we have taken all the masses 
real. In this model, the mass matrix M of the third generation lepton and the vector-like 
lepton is given as following 



jrD-{T,E%)M^ 



(Al) 



and 



^33 ^^34 



II 



E 



(A2) 



V V 

where 777.33 = Vss'^'^^^ ^ "^34 = 2/1'^^ cos /3. Taking ^ I ^34 1; then the biunitary 
matrix to diagonalize Ai'^ are 

LT*j^rj^ri; ^ ^^^^ ^ diag(m^, TJlrl) , (A3) 

where 





( 1 
















V 


1; 





1 




(A4) 



The mass matrix A^^ of the third generation down-quark and vector-like down- type 
quarks is given as following 

CD-{b,Ql,D'k)M' Dl , (A5) 



where 



(m\^ m^4 

77743 "^44 

\ /^"^ 



\ 



(A6) 



where 771^3 = Z/33— p cos ^, 771^4 = y^—j= cos ^, 77143 = Vz' ~j= ^^'^ 



V2 

b I 1^6 



V2 



Taking that ~ \fi^\ 3> |77i4^|, [777441, i"^33h 
Ai^ are 

lh*j^h^b] ^ (^^^^ ^ diag(7n6, 7nbi, 77162) , 



771^4 = y^^-j= COS/ 



7/734!, then the biunitary matrix diagonalize 



(AT) 



15 



where 



and 



1 







1 



\ 



(M^)2+(^Q)2-(m^)2 
1 



(A8) 



1 



(A9) 



The mass matrix of the top quark and vector-hke up-type generations is given as 
following 



CD-{t, Qi U'h) M' 



where 



777.33 ^ 



34 







\Q'hJ 



(AlO) 



777.^3 777^4 jl'^ 



(All) 



where 777*33 = 2/33;^ sin /3, 777I4 = 7/^-^ sin/3, 777^3 = 7/^-^sin/3, 777^4 = y^^-j=sinl3 and 

77i|^ = y—=cosP. Taking that \^^\ ~ \^^\ » |?7743|, I77744I, I777I3I, I77734I, |777|j^|, then the 
v2 

biunitary matrix diagonalize are 



lt*^t^t\ ^ ^^^^ ^[/^ ^ diag(777t, 777ti, 777*2) , 



(A12) 



where 



and 



1 





\ 



1^34 



J 



\ 



(A13) 



(A14) 
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The charged slepton mass-squared matrix of f and the superpartners of the vector- 
hke leptons is given as following 



El* 
El* 



(A15) 



where 



{Ml] 



11 



^ + [-^ - miy] cos2(3 + ml + \m. 



T 12 

34 1 



;i2 



(mo — //tan/3)m^ , 



22 
24 



{Ml) 



(mo - /itan/3)m34 , {MDu = ^i^ml^ , {Ml)2i = {mo - /xtan/3)m^ , 

Me - {'iTT'Z - '^w) cos 2/3 + , {Ml)2z = rrirml^ , 

, (A4^)3i = (mo - //tan/3)m54 , (Al^)32 = m^m54 , 

|^S|2 + M| - (m| - m^) cos2^ + \ml^\'' , {Ml)34 = S^//^ , 

/x^m54 , (yW^)^^ ^ , (MlUs = 



1/^ 



£;i2 



mc 



"z - ^w) cos 2/3 . 
The corresponding unitary scalar matrix is defined as 



(A16) 



(A17) 



The mass-squared matrix A^^ of 6 and the superpartners of the down-type vector- like 
fermions is given as following 

/ b \ 



T-\C* 

^3 

-^4 

Ql 
\Q'hJ 



(A18) 
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where 



b)l3 



_!4jj!!kcos2/3 + m2 + |m^34p , iMl)i2 = (mo - /itan/3)mfe , 



- cos 2/3 + ml + 



= , {Ml)2i = (mo - //tan/3)m6 , {Ml)22 = 3 ■ cuszp -r iiif^ -r \'i 

m6m^4 + m^gm^^ , (^4^)24 = , (^^6)25 = {ruo - /xtan^)m^3 , {Ml)2e = 

( rrir, — II fan R\'m^ , ( A/f?^or. = mim^ . -I- m^.^m^. . 



(A^b)3i = (mo - )Utan^)m|4 , (^4^)32 = m6m^4 + m|3m^4 , 

(Ml):,, = |u^P + M2 -^^^i^cos2/3 + |mtP + |m5J2 



-y^d 1 2 
''■43 1 ) 

//^m^3 , 



(-MD33 
(-MD35 



3 COS 2/9 + |m|4p 



(^46^45 = A*"m44 + ii-'niH , {Jvid^g = l,mo - /icoipjmH , {Jvid^i = m, 
(^4^)52 = (mo - ntanp)m\r^ , {MDbs = {ruo - //tan/5)m^4 , (Alg)54 = 
{Mlh, =|/.«P + M|-^^^i^cos2^+K4P + K3P, (>(^)56 = , 

The corresponding unitary scalar matrix is defined 



(mo - /itan/3)m^4 , {MDse = fi'^m 

, (^,2)44 = I^^P + Mlj, + ^^^i^ cos 2/3 + +|m^p , 
//'°m|4 + , {Ml)46 = (mo - //cot/3)m^ , (TW^jsi = m6m^3 + m^4m^4 , 



'44 



, (-MD34 



, (>lg)63 = /iM4 + Ai^"^/ 



6 12 

■h\ 



+ 



r > (-^6)64 

COS 2/3 



//^m^4 + /i^m^ , 

= (mo — /i cot /3)mff , 
(A19) 



y^>(^y^t ^ diag(M2^ M,\, M,|, M23, m2^, m2^) , 



(A20) 



The mass-squared matrix of f and the superpartners of the up-type vector-hke 
fermions is given as following 

/ i \ 



u 



H 



(A21) 



\Q'hJ 
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where 

{M^^)u =Ml + cos 2/3 + ml + Im^J^ , {^Ml)i2 = (^o - /icot/3)mi , 

{Mfju = {^0 - /xcot/3)m^4 , {Ml)u = l^^rnl^^ , (-^?)i5 = mml^ + ^^44^*4 , 

(A4|)i6 = , (A4|)2i = (mo - Atcot/3)mt , (>(t^)22 = + |(m| - m^) cos 2/3 + + \mi^ 

iMt)23 = mtm|4 + m^3m^4 , (^^4)24 = , (^[4)25 = (^o - iicotp)ml^ , 

(A^f)26 = jw^m^s > (-^Dai = (^0 - cot;5)m|4 , (A^f)32 = Km^4 + ^43^44 , 

{MlU = l/x^P + M2 + |(m| - m2^) cos 2/3 + \mU\^ + |m*4|2 , (TWf)^^ _ -S^/x^ , 

(A^?)35 = (mo - /icot/3)m*4 , (7W?)36 = /u^m*4 + /i^m^ , (7W?)4i = /i^m*4 , 

(>i?)42 = , {Ml),3 = -B^fi'' , (MDu = l/x^P + - |(m| - m^,) cos 2/3 + |m*,|2 , 

(A4f)45 = A*^*m|4 + //<5m^ , (Al|)46 = (mo - //tan/3)m*^ , (TWDsi = mtm^g + m*34mi4 , 

(Mt)52 = (mo - Ijl cot p)m\r^ , (>i|)53 = (mo - //cot^)m^4 , (A4f)54 = /i^mi4 + //'^m^j , 

(A')55 =|//«P + M| + ^^^^^C0S2^+K4P + K3P, (A' )56 = , 

(A^?)6i = , (-M?)62 = yu«m*3 , (yW?)63 = /i'^m*4 + /i^m*^ , (A^?)64 = (mo - /xtan/3)m*^ , 
(A')65 =-5'?/.^, (A')66 = |/x«P + M2^ + |m*,p-^^^^^cos2/3. 

(A22) 

The corresponding unitary scalar matrix is defined as 

V'MjV'^ = diag(M,^ , M,|, M,|, , M,! ) , (A23) 

Appendix B: NEUTRINO MASSES IN OUR MODEL 

AU fourteen types of one loop Fyenman diagrams which can contribute to the neutrino 
mass and mixing in our model are shown in Fig. 1 
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FIG. 7: One-loop contributions to the neutrino masses and mixings in our model. 
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corresponding analytical results are listed below: 
1 

8^ 



k,m 



g^2 / V ^'i3^'JJJ r''ml ml fcl ' fc3 ""Tm^y" "Tmt ■''■^Ti^(R)kJ 
k,m 

k,m 
k,m 

+ A-32Aj-23 RmiKa^br. sin "si(2) cos asi(2)&(m6^, M^^ J] 

fc,m 
k.m 

^5Ia«A- ^ ^E^?3Af3[<2i^;^nW2^6.&K.,M^,(,,J 

k,m 

fc,m 
k.m 



L(R)k ' 



^^■Ia'aq ^ ^E'^^33Ag[i?;.\i^;;.\^;2'^.5^6.K^;>.,M^,^^^^ 

fc,m 

^;^Ia'A- ^ ^E43Af3[^m2^m\^«2^;..K^;>.,^S^^^^^ 

fc,m 



L(R)k' 
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k,m 
k,m 

In which L^'^'*, i?^'^'* are biunitary matrices of mass matrices between (r, 6, t) and the 
vector-hke fermions (see Appendix A), while ^r^,^6m,'^t^ indicate the corresponding 
mass eigenvalues, y^'*^'* are the square mass mixing unitary matrices of their superpat- 
ners, while Mf^^^^^ , M^^^^^^ , M^^^^^^ stand for the corresponding mass eigenvalues. sinQ;si(2), 
cosQ;5i(2) are the unitary matrix elements of s. h{m\,m2) is the loop integral factor: 
6(7711,777.2) = j^-2^j^-2 {ml Inml — m\\!a.m\ — m\ + 777^). The value range of the indices in 
Eq. (4)-(6) is m=l,2, k=l-4 ,while in Eq. (7)-(17), it is m=l,2,3, k=l-6. 



Appendix C: NEUTRINO SPECTRUM-CALCULATING METHOD AND PA- 
RAMETER SETTINGS 



The methods to generate neutrino masses and mixing angles with one-loop trilinear ^ 
couplings actually involves the following three matrices 



7?7i 



a? ah ac 



ah 6^ he 
ac he <? 



\ ( 



\ ( 



9^ gh gl 
gh /i^ hi 
gl hi f 



(CI) 



(P de df 
de ef 

yac be cf j \df ef pj ygi ni r j 

where we name each of the matrices above A1i,2,3 separately. We assume 7?7i > 7772,3, 
77^2 ~ 77^3 and there is no strong hierarchy between a, b, c, d, e, f, g, h, 1. 
A^i has only one eigenvalue after digonalized by an unitary rotation 



where 



and 



X'^MiX = diag(0, 0, Ml) , 



Ml = 777i(a^ + 6^ + c^) , 



(C2) 



(C3) 



^ C2 S2C3 S2S3 ^ 



X 



-S2 C2C3 S2S3 
- S3 C3 



(C4) 



22 



S2 = 



, C3 



Va^ + 62 ' Va2 + 62 + c2' 
If we rotate the sum over 1,2,3 by matrix X, it becomes 

X^{Mi + M2 + M3)X ^ mi(a2 + 6" + c") 



(C5) 



v 



(C6) 



cii €12 ei3 

^21 ^22 ^23 
£31 £32 1 

where e^j are some small values related with 1712/1711, 1713 /mi and the other elements of 1,2,3- 
We can then define another unitary matrix X' to diagonalize the matrix in Eq. (B6) in an 
approximate way: 

X'^X^{Mi + M2 + M3)XX' f«mi(a2 + 62 + c2) diag((5^, Si,, 1) 



where 



— Si Ci 
1 



and 



tan 2^1 



2e 



12 



£22 — fill 

Then from Eq. (B7), we get all three mass eigenvalues 

Ml ~ mi(a^ + 6^ + c^) , M2 ~ MiS, , M3 ~ Mi^g , 
and from Eq. (B4, B8), we get all three mixing angles 

Sl3 = S2S3 = > 
S23 = C2S3/C13 = , 
S12 = (S1C2 + CiS2C3)/ci3 ■ 

The parameter settings we used in table I are given as following 

Set I: 



(C7) 



(C8) 



(C9) 



(CIO) 



(Gil) 



ml^ = 10, Me = Meh = 600GeV, B^/i^ = 400^GeV^; = 170GeV, 
ml^ = 771^3 = = 0, Mq = Md = Mdh = Mqh = 700GeV, 
S^/x^ = B^H^ = SOO^GeV^; w}^ = m^g = = 13GeV, w}^ = 174GeV, 
M[7 = MuH = 700GeV, S^/x^ = bm^GeV'^, tan ^ = 10, A = = 500GeV. 
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Set II: 



ml^ = 10, Me = Meh = lOOOGeV, 5^//^ = GOO^GeV^; = 170GeV, 
m^4 = m^3 = m^4 = lOGeV, Mq = M^^ = Mj^j^ = Mqh = lOOOGeV, 

= = eOO^GeV^; m*4 = = m*^ = 13GeV, m^^ = 174GeV, 

Mu = MuH = lOOOGeV, S^//^ = GOO^GeV^, tan/3 = 10, A = // = 600GeV. 

Appendix D: EXOTIC QUARK AND LEPTON COUPLINGS TO W, Z, hP AND 
DECAY WIDTHS 

The couphngs for the W, Z, boson with leptons in Eq. (12) are 



where = s^, = — is the elements of the rotation matrix related with the real parts 




(Dl) 




r(Ti ^ Zr) 



r(ri ^ WUr) 




(D2) 



r(Ti ^ h°T) 



m. 



■Tl 



(1 + 4o + - 2x^0 - 2x1 - '^xioxiy/^ 




where Xi — mi/mj-^ for i — W, Z, r, h°. 
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The couphngs for the W, Z, hP boson with t, ti, t2 in Eq. (13) are 

9i2LbL = "^(-^31-^11 + -^32-^12) > fi'C^bjj = ;^-R33^13 ) 

_ _9_(Tt jb j_ Tt Tb \ „W _ _9_T3t r>b 

9t2LbiL ~ ^1-^31-^21 "T -^32-^22^ ' ^t^RbiR ^ ^2 33-"'23 ' 



dfiLtL ~ of ^-^23-^13 + Tr~[(2 - q^V^) (-^21-^11 + -^22-^12) 



9hRtR = -^^[■^^w{R2iRii + ^22^12) + (2 - 2Svf)^23^13] 
^ 3^^^ -^33-^13 + 4^['^^ ~ 2 "^vf) (-^31-^11 + -^32-^12)] 

fl'Ifltfl ^ -^^[-^^w{R\lRll + -^32-^12) + (2 - gSw)-^33-^i3] 



9hLtiL ~ 3^^^ -^33-^23 + 4^^^-^ ~ g-^W^) (-^31-^21 + -^32-^22; 
^'IfiiiH = -^^[3"5vf(-R31-R21 + -^32-^22) + (2 - gSw')-^33-^23] ' 

^ (^33-^21-^11 + -^22-^11 + Z/3^ -^21-^12 + ^^^-^22-^12) ~ ~^y^ ^23^13 

9htiR ^ "^(1/33-^11-^21 + Vs^u^ii + 1/3^-^11-^22 + l/?*^ -^12-^22) ~ "^y^^u^ls 

9i2LtR ^ (1/33-^31-^11 + -^32-^11 + 1/3^-^31-^12 + 1/^^-^32-^12) ~ -^y^L33^l3 
9tt2R = ^(1/33^*11^^1 + l/3^^l2i?L + l/3^^*il^k + l/3^^^*i2^k) - ^1/''^*13^^3 
5'f2itlJ{ ~ (1/33-^31 -^21 + l/?-^'32-^21 + y3^il^22 + l/^^-^'L-^D ^ ~^y^ ^33^23 

9hLt2R ^ (1/33-^21-^31 + 1/^-^22-^31 + 1/3^-^21-^32 + 1/3^^-^22-^32) ~ ;^1/'^ -^23-^33 

The decay widths of the hghtest new up-type quark ti are 

r(ti ^ Wb) = ^(1 + + - 2a;2^ - 2a;^ - 2a;2^a;^)^/2 

{(1 + x,^ - 20;^ + (1 - a;D^a;-)[(,,^^,J^ + (,g;,J^] 

+ i2x6^lj;6^^lj;,,^} > 



r(ti ^ zt) = ^(1 + 4 + - 24 - 2x? - 24x^^2 



{(1 + x,^ - 24 + (1 - xlfxf)[{gl,^f + (,f^,J^] 
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where Xi = mi/mt^ for i = W, Zjtjh'^. The heaviest new up-typc quark t2 has six decay 
channels. The decay widths have the similar forms and can be deduced straightforwardly. 
The couphngs for the W, Z, boson with 6, 61, 62 in Eq. (14) are 

S^LtL ^ ^'-■^21-^11 + -^22-^12) > 9^Rtn ^ ~^^23Rl3 > 

r,W — 3 ( rb Tt j_ Tb Tt \ AV _ 9 TDb JDt 

9b2LtL ~ "^/f'^ 31^11 + ^32^12J , 9l^j,tR - ;y|^33«13 > 

r,W — 9 ( j-b Tt , Tb Tt \ „W _ 9 ryb r>t 

9b2LtiL ~ ^^^31^21 ^ ^32^22) > 9b2ntiR- ^-^33^23 > 



^ 4 

9liLbL ^ 3c^ "^23-^13 ~ 4^t^^ ~ g-Sv^) (-^21-^11 + -^22-^12)] 



f^flbfl ^ 4^ ^3^^^- "^21-^11 + -^22-^12) + (2 - gSw')-^23-Rl3] 

„Z _ 9Sw Tb Tb S^r/r, _ Wr6 r6 , r6 rb 

"-^33-^13 7~il^ Q "Sw'J 1-1^31-^11 + -^32-^12jJ 



02 L 



bL 



3cvK 4.cw'" 3 



S'Cflfefl ~ 4^[3'^v^(-^3i-^ii + -^32-^12) + (2 - 3S^)-R33-^13] ' 

^ 4 

S'Cibii ^ 3^-^33-^23 ~ ^^[(2 - gSv^) (-^31-^21 + -^32-^22)] ) (D5) 
9l2RbiR = ^^[3'5vf(-R31-R21 + -R32-R22) + (2 - gSw^)-^^-^^] ' 

= -^(^33^^i^?i + y?Ll,R\, + ya^L^ii?^^ + vTlIA) + ^y^'L^3i?L , 

9tb.R - -^iyi3L'nRli + i/3^'^?2^^i + y?L\,Rl, + 1/3^^^?2^^2) + ^y'^'LlA , 



9hLbR ~ (^33-^31-^11 + -^32-^11 + -^31-^12 + ^^^-^32-^12) + ^^^^ -^33-^^ 

c. 
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9lLb2R ^ ~;^(?/33-^ll-^31 + 1/3^ -^12-^31 + 2/? -^11-^32 + l/s^'^-^ 12-^32) + "^^^ -^13-^33 
^52l6ij{ ^ ~;^(?/33-^31-^21 + -^32-^21 + Z/f -^31-^22 + Z/^'^ -^32-^22) + "^^^ -^33-^23 
^5iL62fl ^ ~;^(?/33-^21-^31 + ?/3^ -^22-^31 + 2/f -^21-^32 + 113^^22^32) + ^^^^ -^23-^33 
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The decay widths of the hghtest new down-type quark bi are 

r(bi ^ Wt) = ^(1 + xt^ + xt- 2xlr - 2x1 - 2xWt)'^' 

{(1 + X? - 2x1 + (1 - ^'tr^w')H,j' + idLj'] 

-.W 



r(bi ^ Zb) = ^(1 + x\ ^x\- 2x\ - 2x1 - 2^14)'^' 

{(1 + xi - 2x1 + (1 - <)Wn9Luf + Ka^J] m 

r(bi ^ h°b) = ^{l + xlo + xt- 2xlo - 2x1 - ^xloxiy/^ 
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{ii+xi-xi.)u:^j^+i9^ij^] 

where Xi = nii/mh^ for index i = W, Zjbjh*^. The heaviest new down-type quark b2 has 
six decay channels, the decay widths have the similar forms and can be deduced straight- 
forwardly. 
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